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ABSTRACT 

Recent results based on the analysis of radio galaxies and their hot X-ray emitting atmospheres 
suggest that non-radiating particles dominate the energy budget in the lobes of FRI radio galaxies, 
in some cases by a factor of more than 1000, while radiating particles dominate the energy budget 
in FRII radio galaxy lobes. This implies a significant difference in the radiative efficiency of the 
two morphological classes. To test this hypothesis, we have measured the kinetic energy flux for 
a sample of 3C FRII radio sources using a new method based on the observed parameters of the 
jet terminal hotspots, and compared the resulting Qjot — ^radio relation to that obtained for FRI 
radio galaxies based on X-ray cavity measurements. Contrary to expectations, we find approximate 
agreement between the Qjet — -^radio relations determined separately for FRI and FRII radio galaxies. 
This result is ostensibly difficult to reconcile with the emerging scenario in which the lobes of FRI and 
FRII radio galaxies have vastly different energy budgets. However, a combination of lower density 
environment, spectral ageing and strong shocks driven by powerful FRII radio galaxies may reduce 
the radiative efficiency of these objects relative to FRIs and couteract, to some extent, the higher 
radiative efficiency expected to arise due to the lower fraction of energy in non-radiating particles. An 
unexpected corollary is that extrapolating the Qjet — ^radio relation determined for low power FRI 
radio galaxies provides a reasonable approximation for high power sources, despite their apparently 
different lobe compositions. 
Subject headings: galaxies: active - galaxies 



1. INTRODUCTION 

Measuring the kinetic power of extragalactic jets 
has application in two important areas of astrophysics: 
(1) determining the Active Galactic Nucleus (AGN) 
kinetic luminosity function and its evolution — an 
impor tant factor in the s tudy of radio galaxy feed- 
back (Croton et alj |2006|; iShabala fc Alexandeij 120091 : 
iFanidakis ct al. 20 111): and (2) assessing the contribution 
of bl ack hole spi n in the production of extragalactic jets 
(e. g. iDaivl 120091 and references therein). Furthermore, 
accurate estimates of jet kinetic energy flux can provide 
an important constraint in assessing X-ray emission mod- 
els of kpc-scale quasar jets (Godfrey ct al. 2012). Accu- 
rate measurement of jet power in radio galaxies is key to 
qua ntifying their effect on ga laxy evolution. For exam- 
ple, iRawlings fc JarvisI ()2004[ ) and IShabala et all ()2011f ) 
showed that bow shocks driven by the high power radio 
sources can suppress star formation in not only the AGN 
host, but also in other nearby galaxies. More impor- 
tantly for the current work, the ratio of radio luminosity 
(^radio) to jct powcr (Qjet), slso known as the radiative 
efficiency, is sensitive to the division of energy between 
radiating and non-radiating particle populations, and 
therefore may be used to investigate the lobe energetics. 
However, measuring the kinetic power of radio galaxies 
has proven to be a very difficult problem. The lack of re- 
liable empirical methods to measure the kinetic power of 
AGN jets has resulted in the widespread use of a model- 



dependent predictor of jet power derived bv lWillott et al.l 
(|1999| ). based on synchrotron minimum energy calcula- 
tions in combinat ion with th e self-similar model of radi o 
galaxy evolut i on (lFallelll991l : iKaiser fc Alexanded [19971 ) . 
IWillott et all 119991) obtain an expression for the jet 
power Qw ("W" for Willott) in terms of the 151 MHz 
radio luminosity 
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where Qw is time averaged kinetic power of a source with 
radio luminosity L151 = Fi^iD\, and / is a parameter 
accounting for systematic error in the model assump- 
tions. These model assumptions include, among other 
things, the fraction of energy in non-radiating particles, 
the low frequency cutoff in the synchrotron spectrum, 
and departures from minimum energy. It is argued by 
IWillott et all (|1999D that 1 < / < 20, implying a system- 
atic uncertainty of 2 orders of magnitude in jet power for 
a given radio luminosity, owing to the /^^^ dependence. 
The Willott et al. jet power relation is widely used to 
estimate the mechanical output from AGN based on a 
single low frequency luminosity measurement, assuming 
that the value of / is constant (typically of order 10 - 
^0) across the entire population of radio galaxies (c. g. 
Hardcastlc ct al. 2007; Martincz-Sansigrc & Rawlingi 
20111: iFernandes et al.l 1201 It iCattaneo fc Best) l2009). 
The value of / is often calibrated against FRI radio 
galaxies for which the jet power has be en determined 
based on the observed X-ray cavities (e. g. lRaffertv et al.l 
120061) . However, this procedure of calibrating the value 
of / for FRII radio galaxies based on measurements 
of FRI radio galaxies may not be appropriate because 
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of their vastly different energy budgets: non-radiating 
particles are thought to dominate the energy budget of 
FRI radio lobes by a fac t or of ^ 100 in s ome cases 
(jCrostonet all [2001 [200l iBirzan et al.1l2008f ). while in 
FRII radio galaxies, radiating parti cles are thought to 
dominate the ener gy budget (CrostonjrtjilJ [2004, 2005; 
iBelsole et al.ll200"7l ). indicating that significantly different 
values of / should apply to the different morphological 
classes. Indeed, it appears that vastly different values 
of / apply to FRI radio galaxies with different evolu- 
tionary histories (Cavagnolo ct al. 2010). Moreover, the 
self-similar model of radio source evolution, on which 
Equation [T] is based, does not strictly apply to FRI radio 
galaxies. 

Not only is the normalisation of the iWillott et al.l 
(Il999f) relatio n high ly uncertain, so too is the exponent. 
IWillott et all {1999)' use the Oil narrow line luminosity 
(Lou) as a proxy for jet power, and argue that Equation 
[jis valid because the power law exponent matches that of 
the Lou — Li^i correlation that they find in their sample. 
However, the slope of the Lqii^Li^i correlation is highly 
uncertain, and strongly depends on the sample involved: 
iHardcastle et al.l p009.) find Lqii oc Llf^^"-^, while 
iFernandes etall (|20Tll) find Lq// oc L°fs^^-\ Moreover, 
Lqii is expected to exhibit a fairly weak dependence on 
accretion power, and the relationship is likely to be non- 
linea r i n general ([Tad huntcr ct al 1998: Hardcastl e et al.l 
I2009| ). iHardcastle'et al. (2009) find that the correlation 
between Lqii and accretion related X-ray emission is 
not significant after the com mon correlation with red- 
shift is accounted for, however. iShabala et al.l(|2012[ ) find 
a correlation between jet power and narrow line luminos- 
ity for a sample of flat -spectrum radio quasars. Finally, 
lO'Sullivan et al.l (|2011[ ) showed that when the minimum 
energy calculations are performed correctly, the expo- 
nent in Equation [T] depends on the spectral index of the 
source, which in turn is a function of source age (Shabala 
& Godfrey 2012, hereafter Paper 2; Blundeh et al. 1999). 
In the present work, we seek to address these outstand- 
ing questions regarding the Qjct — ^radio relation in FRII 
radio galaxies using an empirical means to determine the 
jet power. 

Jet power in FRII radio galaxies and quasars can be 
estimated directly from measurements of the hotspot size 
and equipartition magnetic field strength, along with a 
number of reasonable assumptions regarding the hotspot 
plasma. In this paper we provide a derivation and anal- 
ysis of this technique, and investigate the relationship 
between jet kinetic power and monochromatic radio lu- 
minosity for a sample of high power FRII radio galaxies. 
We compare the resulting Qjct — iradio relation to the 
equivalent relation obtained for FRI radio galaxies based 
on X-ray cavity measurements, as well as that derived 
by Willott et al. (Equation [T]). 

In Section 2 we provide a derivation of the jet power 
measurement technique based on the parameters of the 
terminal hotspots. In Section 3 we describe the sam- 
ple selection and analysis. In Section 4 we present our 
results, including a comparison of the Qjet — ^radio rela- 
tions for FRI and FRII radio galaxies. In Section 5 we 
discuss the implications of our findings and in Section 6 
we present our conclusions. 

Throughout this paper, the spectral index (a) is 



defined such that Si, oc and we adopt the 

following values for cosmological parameters: Hq — 
71 km Mpc-\ flm = 0.27, Ha = 0.73. Note that, 
following convention, we define the radio luminosity as 
as Li, = Si,D\, where 5*^ is the fiux density, and Dl is 
the luminosity distance. 

2. JET POWER IN FRII RADIO GALAXIES BASED ON 
MEASUREMENTS OF THE TERMINAL HOTSPOTS 

Consider a uniform jet of area A, particle energy den- 
sity e, pressure p, mass density p, relativistic enthalpy 
w = e -\- p + p(? ^ magnetic field components perpendicu- 
lar and parallel to the flow direction i?^ and i?| | , speed 
/3c and corresponding bulk Lorentz factor F. The fiux of 
energy {Fp, ) and rn omentum (Fm) along the jet are (e. g. 
IDouble et aD[200l : 
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In a highly relativistic jet with F ^ 1, the energy flux 
-pEjot is simply related to the momentum flux Fjvijot via 
(see the Appendix) 

^E,jot«C X FMJet- (4) 

Conservation of momentum between the jet and hotspot 
then implies 
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where -pM.hs is the momentum flux in the hotspot. 
This equality holds regardless of assumptions about the 
jet characteristics such as its composition or the ratio 
of magnetic to particle energy densities in the jet or 
hotspot. Therefore, if the general principle of conserva- 
tion of momentum applies between the jet and hotspot, 
we can estimate the jet kinetic luminosity simply by cal- 
culating the momentum flux in the hotspot. To do so, we 
assume that the lepton population in the hotspot is ultra- 
relativistic {e^± = 3pe± ) and that the proton population 
in the hotspot, if it exists, is at best mildly relativistic, 
and can be approximated as a thermal gas (e^ = (3/2)pp) 
so that the hotspot enthalpy density, w, is parametrized 
as follows 
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We further assume that in the hotspot, the magnetic fleld 
is aligned perpendicular to the jet direction, i. e. B± = B 
and = 0, consistent with radio polarisation maps of 
hotspots. The hotspot momentum flux is then: 

+ AceB [l + 2r2/32] 

Let Bcq be the equipartition magnetic fiel d strength, ca l- 
culated using standard expressions (e. g. IWorralll[200l . 
assuming negligible energy density in non-radiating par- 
ticles. Without loss of generality, we can write 
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and then express the jet energy flux in terms of 
Qhs = a c X g 



(8) 
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2.1. Empirical determination of g, the normalization 

factor 

In this section, we empirically determine the value of g 
in equation |8] by applying the hotspot method to sources 
with independent jet power measurements. 

The prototypical FRII radio galaxy Cygnus A is the 
prime candidate for calibrating the normalisation fac- 
tor, since its hotspot parameters are well determined and 
the jet power has be en independently m easured using a 
variety of methods. IWilson et al.l (j2006D estimated the 
jet power of Cygnus A to be Q > 1.2 x 10''^ erg s~^ 
based on analysis of the cocoon dynamics d etermined 
using Chandra X-ray imaging spectroscopy. Ilto et all 
(|2008f ) obtain a similar estimate of jet power (Q = 
0.4 — 2.6 X 10^^ ergs s~^), based on dynamical modelling 
of the source. An independent estimate of jet power in 
Cygnus A comes from Lobanov (19_9j) who shows that 
frequency dependent shifts of the radio core enable a 
deter mination of the jet power (see also iShabala et al.l 
I2012D , and when applied to the case of Cygnus A, gives 
Q « 6 X 10*^ erg s'^ . iRaff ertv etall (|2006f ) estimate the 
jet power of Cygnus A from the X-ray cavity, and find 
Q « 1.3 X 10**^ erg s~^, but more recent analysis incorpo- 
rating the energy in shocks suggests Q > 5 x 10"*^ erg s~^ 
(P. Nulsen, Priv. Comm. 2012). We estimate the jet 
power of Cygnus A using Equation [5] along with the pa- 
rameters fo£_thfi_tenmi]B^ A and D) 
given m iWright fc Birkinshawl (|2004[ ). The derived jet 
power is (5 = g X 4 X 10''^ ergs s~^ , which implies g Ki 1 — 2 
in this source. 

We can apply a similar analysis to the FRII radio 
galaxy 3C 401, as it has an independent jet power mea- 
surement from analysis of the associated X-ray cavity: 
[Raffertv et al, (2006) estimate the jet power of 3C401 
to be Q « 7 X lO''^ erg s~^. This source does not ap- 
pear in our sample because the beam size of the high- 
est resolution map available to us is larger than 1% 
of the source linear extent, and therefore, the hotspots 
may not be adequately resolved (see Section [3]) . The 
hotspots in this source are faint relative to the sur- 
rounding lobe emission, making the determination of 
hotspot parameters difficult. Despite these problems, 
based on th e hotsp ot sizes and flux density given by 
iMuUin et al.1 ((200.81), we obtain a hotspot jet power es- 
timate oi Q — g X 3.2 x lO"*** ergs s~^, which implies 
(7-2. 

Finally, we note that lDalv et al.l (|2012D estimate the jet 
power for a sample of 31 high power FRII radio galaxies, 
using the expression Q = ApV/r, where p is the lobe 



pressure calculated using minimum energy arguments, V 
is the lobe volume assuming cylindric al symmetry , and r 
is the spectral age of the source. The lDalv et all ([20l2h 
sample includes four sources from our sample (3C 55, 3C 
244.1, 3C 289 and 3C 337). For these four sources, the 
av erage value oi g re quired to match the jet power values 
of iDalv et all (|20T1 is g = 2 ± 1. 

The low value derived for the normalisation factor {g — 
2) indicates that the hotspot plasma is close to equipar- 
^+°') tition conditions, consistent with the results of inverse 
Compton niodc]ling in hotspots (Hardcastlc et al. 200.^ 
and lobes (|Croston et al]i2005l) of FRII radio galaxies. 

(9) 

2.2. Expectations for g 

Here we argue that the value for the normalisation fac- 
tor derived from observations in the previous section is 
consistent with expectations, by combining constraints 
on the various parameters involved in Equation |9l 

2 

We first consider the ratio We can constrain this 

ratio, assuming at most one proton per radiating lepton, 
via an estimate of the mean electron Lorentz factor (7) 
in the hotspot, since 
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For a power-law electron energy distribution of the form 
iV(7) = fce7^", the mean electron Lorentz factor (7) w 



7min (fr^ j for a > 2 and (7) = 7niin In [j^j for a = 2. 
In each case where flattening of the hotspot radio spec- 
trum has been directly observed, estimates of 7min are 
in th e order of several hu ndred: PKS 1421-490, 7,„i„ — 
650 (IGodfrev et al.l[ 2?)09'): Cygnus A, 7min ^ 3 00 - 40 
(iCarilli et al JllMlt tazio et al. 2006; Hardcas tl3l2001bD: 
3C295 7min ~ 800 (|Harris et al. 2000; Hardcastld 
2001bl): 3C123. -/ ^in ~ 1000 j Hardcastle e t al. 2001 J 
Hardcastld r2001bf) . In additionTEardcastld (|l2001b ) in- 
ferred a cutoff Lorentz factor 'y^i^ ~ 500 in 3C 196 by 
synchrotron self-Compton modelling of the hotspot spec- 
tral energy distribution. All of the above listed 7niin es- 
timates appear to be distributed around a value of order 
7min ~ 500, and this value of the minimum Lorentz fac- 
tor may arise naturally throu gh the dissipation of bulk 
kinetic energy in the hotspots (IGodfrev et al.ll2009f ). We 
therefore assume 7min ^ 500, and with typical elec- 
tron energy index a ^ 2.2 — 2.6, Equation (TOj implies 

pcVee± < 1. , , 

We now consider the ratio ep/e^±. iHardcastle et al.l 
(|200l argue, based on synchrotron self Compton mod- 
elling of hotspot X-ray emission, that an energetically 
dominant proton population is disfavoured, such that 
£p/ce± ^ 1- The lobes of FRII radio galaxies are in- 
flated by backflow of hotspot plasma. It has been shown 
that ec/ep± < 1 in the lobes of FRII radio galaxies 
(|Croston et al.ll200l I2005L IBelsole et al.]l2007l) . further 
suggesting that the ratio ep/cgi < 1 in hotspots. 
We next consider the ratio IHardcastle et al.l 

(pool argue, on the basis of synchrotron self Compton 
modelling of X-ray emission in a large sample of FRII 
radio galaxy hotspots, that magnetic field strengths are 
close to the equipartition estimates; that is 7^ ~ 1. 
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Fig. 1. — The normalization factor g (see Equation [Sj as a func- 
tion of B/Bcq for a range of assumed states of the hotspot plasma. 
Without loss of generality, we have defined Bcq to be the equiparti- 
tion magnetic field strength calculated assuming negligible energy 
density in non-radiating particles. The solid curve corresponds to a 
case in which the jet is purely leptonic and the hotspot plasma has 
a non-relativistic velocity. These assumptions correspond to the 
minimum possible jet power for a given hotspot size and luminos- 
ity, with g = 1.06, at B/Bcq ~ 2/3. The dashed curve corresponds 
to a hotspot in which the radiating and non-radiating particle pop- 
ulations have equal energy density, the rest mass energy density 
is equal to that of the radiating particle energy density, and the 
hotspot maintains a mildly relativistic velocity of 0.33 c. 

Finally, we consider the bulk velocity of the hotspot 
plasma, p. The post-shock velocity following a normal 
shock in an unmag netised plasma with relativist ic equa- 
tion of state is |c. iDennett-Thorpe et al.l ()1997| ) discov- 
ered a hotspot spectral index asymmetry in the sense 
that hotspots fed by the approaching jet have a flat- 
ter spectral index than hotspots fed by the receding jet. 
This observation can be explained if the hotspots are 
moderately Doppler beamed and have curved spectra, 
with higher frequencies corresponding to a steeper spec- 
trum. The observed hotspot spectral index asymmetry 
requires only mod erate velocity in t he hotspot regions, 
with P > 0.3 (Den nett-Thorpe et al.f 1997) . 

Given the constraints on hotspot parameters discussed 
above, we plot the function g to determine physically re- 
alistic values (Figure [TJ. In this Figure we plot g as a 

function of B/B^^ for 1/3 < < 2, < f| < 1, 

Cp/cei < 10, < /? < 0.33, and a = 0.5. It is clear that 
the emprically derived value of g 1 — 2 is consistent 
with expectations. We note that 5^3 for a hotspot 
near minimum energy with a proton/electron composi- 
tion, /3 ^ 1/3 and Ce ~ ~ pc^, while 5 w 1.5 for 
a purely leptonic hotspot near equipartition conditions 
and /3 ~ 1/3. 

3. SAMPLE SELECTION AND DATA ANALYSIS 

We select a subset of the complete, flux limited sam- 
ple of 3CRR FRII radio sources with redshifts z < 1 
(jMullin et all [2008). In order to limit the effects of ob- 
serving resolution, we excluded objects for which there 
are fewer than 100 beam widths across the source in the 
highest resolution maps available to us. Hotspot sizes are 
typic ally 1 percent of the source size (jHardcastle et al.l 
I1998D . and therefore a lower limit of 100 beams along 
the source ensures that the beam size is smaller than, 
or approximately equal to the hotspot size. A total of 



62 sources met this criteria. We further restricted our 
sample to include only those sources for which a clear jet 
termination could be determined for at least one side of 
the source, resulting in a sample size of 30. Using the ra- 
dio maps available from the 3CRR FRII online database, 
we independently measured the hotspot sizes and flux 
densities in a consistent manner: we fit a 2D-Gaussian 
function to each hotspot using the CASA task IMFIT, 
but only included in the fit those pixels with a value 
greater than half the hotspot peak value. This typically 
resulted in 50 — 80 pixels being used in the fit. We do 
not attempt to model or subtract the background con- 
tributed by the surrounding lobe emission, but we have 
mitigated the effect of the lobe emission by restricting the 
model-fitting routine to include only the bright part of 
the hotspot. Also, our sample was chosen to include only 
those hotspots that were well resolved from the surround- 
ing lobe emission, and we expect that for the majority 
of hotspots in our sample, the lobe emission has a neg- 
ligible impact on our derived hotspot size/fiux density. 
For each source, we used the highest resolution FITS im- 
age available from the online database. The hotspot jet 
powers derived in this way are on average 50% greater 
than the values determined u sing the measured h otspot 
size and flux density given bv lMuUin et"all ([20081 ). In at 
least one instance, the fitting routine gave unreasonable 
results, and we then estimated the hotspot region from 
the map, and extracted the flux using an elliptical aper- 
ture of the same dimensions as the hotspot dimensions 
estimated from cross-sectional profiles of the hotspot. 

As with iMullin et al.l (|2008[ ). we have not quoted er- 
rors for hotspot flux density or size since the dominant 
source of error comes from the ambiguity in defining the 
hotspot region, and is therefore, to some extent, sub- 
jective. However, we have minimised the ambiguity by 
restricting our sample to include only those sources for 
which the beam width is less than, or com parable to the 
expected hotspot size (jMullin et al.ll2008l) . and therefore 
expect our derived hotspot powers to be reliable esti- 
mates, free of large systematic bias. 
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Fig. 2. — Jet power, determined using the hotspot method, as a 
function of 151 MHz monochromatic radio luminosity. The solid 
black line is the best fit relation to our data, which are plotted 
without multiplication by the normalization factor (effectively as- 
suming g = 1). 
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4. RESULTS 

For each hotspot, we compute the equipartition mag- 
netic fie l d str ength Bcq using standard expressions 
(IWorralll 120091) assuming negligible energy density in 
non-radiating particles. We assume a power-law elec- 
tron energy distribution of the form N{'-f) = fce7~° with 
a = 2.2 between a minim um Lorentz factor 7i„in = 500 
(see IGodfrev et all 12009 1) and maximum Lorentz factor 
7max = 10^. For each hotspot, we then calculate the jet 
power using Equation [8] with g = 2. We sum the power 
derived from the hotspots on each side of the source to 
obtain the total source power. In Section [5.11 we discuss 
the reliability of jet power measurements from individual 
hotspots, and find that typically the jet power estimated 
for a pair of hotspots in the same source agree to within 
a factor of 2. In a few cases, only one hotspot could be 
used for a reliable jet power estimate, and in that case, 
we multiplied the derived value by a factor of 2, to ac- 
count for the power in the oppositely directed jet. The 
sample and derived jet powers are given in Table [T] We 
perform least-squares minimisation in log space to fit a 
power law of the form Q = ALfg^ to the data, and find 
a best-fit relation 
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Calculating the hotspot jet powers usi ng the measu red 
hotspot parameters tabulated bv iMullin et al.l (|2008| ) re- 
sults in jet powers that are on average a factor of 1.5 
lower. Assuming 7min = 10 rather than 7min = 500 in- 
creases the normalisation of the best-fit curves by a fac- 
tor of approximately 1.5. The assumed hotspot spectral 
index affects the slope of the derived best-fit relation, 
however the effect is comparable to the statistical uncer- 
tainty in the slope due to the scatter in the correlation. 
For example, assuming that the hotspot spectral index is 
a = 0.8 instead of a = 0.6 results in a marginally flatter 
best-fit relation, with an exponent of 0.62 ± 0.05. 

4.1. Comparison with the Qjct — iisi relation for FRI 
radio galaxies 

A number of authors have investigated the rela- 
tionship between AGN jet power and radio power in 
nearby low-luminosity sources using X-ray cavity mea- 
surements; tlie_jnas^_wWfily discussed rec ent work be- 
ing that of iCavagnolo et aL (|2010D (see also lBirzan et al.l 
l2008HO'Sullivan et al. II 20111) . Studies of jet kinetic power 
based on X-ray cavity measurements are inherently lim- 
ited to relatively nearby, low power objects, typically of 
FRI morphology. We seek to compare our results to the 
Qjet — ii5i relation obtained for FRI radio galaxies. To 
do so, we ha ve used the X-ray cavit y jet power estimates 
compiled bv lCavagnolo et aLl ()2010[ ) along with 151 MHz 
radio luminosities extrapolated from low frequency mea- 
surements, assuming a spectral index a — 0.8. The low- 
frequency luminosities for the Cavagnolo sample lie be- 
tween 200 - 400 MHz, and are mostly at 327MHz. This 
is sufficiently close to 151 MHz for the assumed value of 
spectral index to be largely insignificant: a departure in 
the assumed spectral index of Aa = 0.3 results in only 
a ~ 20% error in L151. We have converted the best-fit 



relation of ICavagnolo et ail (|2010t ) to one involving L151: 
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The exponent in the above relation for FRI radio galax- 
ies is in excellent agreement with the exponent deter- 
mined for our sample of high power FRII sources (Equa- 
tion [TT]). We pl ot Equation IT^ i n Figu re |2| along with 
the data for the ICavagnolo et al.l ()2010[ ) sample of low- 
luminosity radio galaxies (red points), and the data for 
our sample of FRII radio galaxies (green squares) using 
the assumption g = 2 (see Section I^TTj) . We find no evi- 
dence for a significant offset between the FRI and FRII 
Qjet — iradio relations: the normalisations formally agree 
if 51 > 2, which is entirely consistent with o ur analy- 
sis m Sections O and O iDalv et al.l ([2012^ estimate 
the jet power for a sample of 31 high luminosity FRII 
radio galaxies using the expression Q = ApV/r, where 
p is the lobe pressure calculated using minimum energy 
arguments, V is the lobe volume assuming cylindrical 
symmetry, and r is the spectral age of the source. In 
agreement with the results presented here, Daly et al. 
also found that the Qjct — ^radio relation for luminous 
FRII radio galaxies is in broad agreement with an ex - 
trapolation of the one given bv ICavagnolo et al.l (poTol ) 
for FRI radio galaxies. 



• • FRI Jet Power Estimates from X-ray Cavities. 

■ ■ FRII Jet Power Measurements from Hotspots ((/ = 2). 

- - Willott relation with /=20 (upper) and / = ! (lower). 

— Extrapolation of FRI Q-L^^^ relation. 




. Q„ (f=20) 



Q„- (f=l) 



24 26 28 30 

Li5i [ W Hz 1 Sr i ] 

Fig. 3. — Comparison of the CJjct — iisi relations for FRI and 
FRII ra dio galaxies. Here we plot the data and best fit relation 
from ICavagnolo et al. (2010) (red points and blue solid line). The 
shaded area illustrates uncertainty in the normalisation of the FRI 
best fit relation. We also plot the model of Willott et al. ( 199J|) 
with / = 20 (upper-most black dashed line) and / = 1 (lower-most 
black dashed line). We plot the FRII jet power measurements 
(green squares) which have b een deri ved using the hotspot method 
assuming g = 2 (see Sections 12. Il and l2. 21 1. Note that the minimum 
allowed value is g = 1.06. The black cross marks the location of 
Cygnus A and is clearly an outlier when compared to our sample 
of FRII radio galaxies. This is due to the high density environ- 
ment into which Cygnus A expands, resulting in "environmental 
boostin g" o f its radio luminosity (iBarthel fc Arnaudiil996i , see also 
Section [5^. 

Our result, the broad agreement between the FRI and 
FRII Qjct — ^radio relations, appears at odds with the 
emerging scenario in which the fraction of energy in 
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non-radiating particles differs greatly between these two 
classes of radio galaxy. However, as we discuss in Section 
15.21 differences in the age and environment for the two 
samples used in this study, as well as a possible differ- 
ence in the fraction of energy associated with shocks, will 
counteract the offset between the Qjet — ^151 relations ex- 
pected to arise due to the differing energy budgets of the 
radio lobes. 




0.3 0.41 0.55 0.74 1.0 1.35 1.82 



Ratio of hotspot derived jet power (with g=2) to that predicted 
from the Q.^^ - L^^^ scaling relation for FRI radio galaxies. 

Fig. 4. — Histogram of the ratio between hotspot j et p ower (with 
g = 2) and the jet power calculated from Equation 1121 the Qjet — 
L151 scaling relation for FRI radio galaxies. The mean of this 
distribution (0.8) is illustrated by the dashed line. It is clear that 
given 3 Si 2 as derived in section 12.11 there is no evidence for a 
substantial offset between the FRI and FRII Qjct — ^radio relations. 



5. DISCUSSION 

5.1. Hotspots as calorimeters 

Hotspots of FRII radio galaxies are thought to be vari- 
able on short timescales (jLai ng 1989; iSaxton et al.ll2002l 
|201C1| ). and as such, caution must be exercised when in- 
terpreting the derived jet power for individual objects. 
However, provided that the general principle of conser- 
vation of momentum applies between jet and hotspot, 
on a population basis we expect this method to be a re- 
liable estimator of jet power, and in particular, may be 
used to investigate the Qjct — iradio relation at high ra- 
dio luminosities. More than half of the sources in our 
sample have two hotspots, one at each end of the source, 
that enable jet power estimates. We can test the reli- 
ability of the hotspot jet power method by calculating, 
for each source, the ratio of jet power determined for the 
two hotspots. Figure[5]is a histogram showing the distri- 
bution of this Qhs ratio. More than half the sample have 
Qjot estimates from both hotspots that agree to within a 
factor of two. The largest discrepancy between hotspot 
measurements is approximately a factor of 5. 

5.2. Predicted offset between the Qj^t ~ -^ladio relations 
for FRI and FRII radio galaxies 

lO'SuUivan et all (|201lD revised the analysis of 
iWillott et al.l ( 19990 to account for a different minimum 
energy formalism. In particular, these authors pointed 
out that a large fraction of the synchrotron-emitting elec- 




Fig. 5. — Histogram of the ratio of jet power derived from the two 
hotspots at either end of the source. The median of the distribution 
is 2.0 and standard deviation is 1.4. 

tron population may radiate below the frequency cut- 
offs assumed bv IWillott et al.l ()1999l ). Because of this, it 
makes more sense to recast the minimum energy argu- 
ment in terms of cutoff Lorentz factors of the electron 
energy distribution. An additional advantage of this ap- 
proach is that it allows the Qjct — ^radio relation to be 
expressed in terms of the lobe spect r al ind ex. Below, we 
briefly recount the lO'SuUivan et all ()2011| ) analysis. 

Let Uniin be the minimum energy of a synchrotron emit- 
ting source of a given volume and luminosity, k the ratio 
of energy in non-radiating particles to the energy in radi- 
ating particles, and i?nio the mi nimum energ y mag netic 
field strength. In the model of iWillott et al.l (|1999[ ) the 

jet power is shown to be Qjct cx: "u^n (see their Equation 
11). The minimum energy Umin is related to the min- 
imum energy magneti c field stregnth via Umin oc B^^. 
IWorrall &: Birkinshawl (|2006l ) provide an expression for 
the minimum energy magnetic field strength in terms of 
the observed source luminosity and volume, as well as the 
high and low energy cutoff in the electron energy distri- 



bution, which has the form B„ 
Combining these expressions 



oc {1 + A:) 3+= L 



l/(3+a) 



151 



Qjetoc(l + fc)^L?/f+") 



(13) 



In Paper 2, we extend this analysis to show that the 
Qjet — iiadio relation is sensitive to a variety of radio 
source parameters, including lobe size and age. It is clear 
from Equation 1 1 3! however . that the Qjet — iradio relation 
will be altered if either the the fraction of non-radiating 
particles k, or the spectral index a change (as expected 
from straightforward synchrotron ageing arguments, e.g. 
Alexander & Leahy 1987). 

We can estimate the expected offset between the Qjot — 
.^radio relations for FRI and FRII radio galaxies due to 
the presence of non-ra diating particles using Equation 
[T51 iBirzan et al.l ()2008l) estimate k for a sample of mostly 
FRI radio galaxies by equating the internal lobe pressure 
with the external pressure, and requiring that the mag- 
netic field strength be in energy equipartition with the 
particles. They find k lies in the range ~ 1 — 4000, with 
a med ian value of k w 180. Similarly, iCroston et al.l 
(|2008f) find that the radiating material in a sample of 
FRI radio lobes is significantly under-pressured relative 
to the external environment, and infer the presence of 
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TABLE 1 

HOTSPOT Properties and Derived Jet Power 



Source L151 Hotspot Beg Dhs'" Q.jct'' (g=2) Qtotal'' (g=2) 

xlO^^ North South North South North South 

W Hz-i Sr-i ^jG ^lG kpc kpc xlQ-'^ergs-i xlO'^^crgs-i xlQ-*^ erg s"! 



3C22 


46.5 


680 


320 


0.93 


3.7 


7.2 


24 


31 


3C33.1 


1.07 


24 


10 


19 


19 


3.5 


0.65 


4.1 


3C46 


7.75 


99 


27 


2.3 


12 


0.94 


1.9 


2.8 


3C55 


55.6 


290 


160 


1.7 


4.4 


4.5 


9.5 


14 


3C98 


0.093 


25 


30 


3.3 


2.3 


0.12 


0.09 


0.21 


3C109 


6.15 


93 


191 


2.5 


2.3 


0.96 


3.6 


4.6 


3C132 


1.65 


320 


105 


0.59 


1.6 


0.65 


0.49 


1.1 


3C184.1 


0.43 


33 


68 


5.0 


2.2 


0.49 


0.43 


0.91 


3C228 


27.2 


280 


260 


1.9 


2.6 


4.9 


8.2 


13 


3C234 


2.88 


170 


120 


0.81 


2.2 


0.36 


1.3 


1.7 


3C244.1 


12.5 


280 


145 


1.8 


2.2 


4.7 


1.9 


6.6 


3C284 


1.92 


30 


81 


6.0 


2.7 


0.60 


0.89 


1.5 


3C337 


16.3 


190 


275 


2.3 


3.1 


3.5 


13 


17 


3C340 


23.7 


91 


250 


9.3 


3.1 


12.9 


10 


24 


3C382 


0.14 


41 


25 


3.8 


4.5 


0.43 


0.23 


0.66 


3C33 


0.43 




160 




0.95 




0.41 


0.82 


3C42 


5.87 




100 




4.6 




3.8 


7.5 


3C223 


0.67 




17 




10.5 




0.58 


1.2 


3C226 


45.3 




470 




1.2 




5.9 


12 


3C263 


25.3 




500 




1.7 




13 


27 


3C277.2 


34.0 




360 




2.2 




12 


23 


3C289 


50.9 




950 




0.67 




7.2 


14 


3C292 


20.3 




130 




5.4 




9.0 


18 


3C300 


3.82 




93 




3.2 




1.6 


3.1 


3C319 


1.57 


42 




4.2 




0.56 




1.1 


3C349 


1.49 




130 




2.7 




2.1 


4.1 


3C381 


1.11 


170 




1.15 




0.70 




1.4 


3C452 


0.84 


27 


69 


8.0 


1.7 


0.87 


0.26 


1.1 


3C321 


0.29 


34 


84 


3.2 


3.1 


0.21 


1.2 


1.4 



^ Hotspot diameter, normal to jet direction. 

^ Jet power derived from hotspot parameters using cquation[S] assuming the normalisation factor p — 2, as derived in Section 

^ Total source power (sum of north and south jet power, or twice the measured jet power if only one hotspot could be used), 
again, assuming g — 2. 



non-radiating particle population that st rongly domi- 
nates the lobe energy budget. In contrast. iBelsole et al.l 
((21307) find that the equipartition lobe pressures of FRII 
radio galaxies are typically close to pressu re equilibrium 
with t he external medium. Furthermore, ICroston et alj 
(|2005l ) find that the magnetic field strength in the lobes 
of FRII radio galaxies, determined via synchrotron and 
inverse Compton modeling of the radio to X-ray spec- 
tra, is close to the equipartition value determined from 
the radio data alone assuming negligible energy in non- 
radiating particles. They conclude that FRII lobes are 
unlikely to contain an energetically dominant population 
of non-radiating particles, because that would require the 
magnetic field energy to be matched to the energy of 
just the relativistic electron population rather than the 
energy of the entire particle population. The inferred 
difference in lobe energy budget is not altogether sur- 
prising, since FRII lobes are surrounded by bow shocks, 
and the jets are surrounded by cocoon plasma, making 
significant entrainment of ambient material rather more 
difficult than in FRI lobes, for which this is not the case. 

If k ^ 200 in FR I radio galaxy lobes (the median value 
from lBirzan et ai1 [20081 and A: ~ 1 in FRII radio galaxy 
lobes, then naively applying Equation [T3] to both mor- 
phological classes with a — 0.8, we would expect the 
normalisation of the Qjct — ^radio relation for FRI radio 
galaxies to be greater than that of FRIIs by a factor of 



> 50. This is clearly not observed (see Figures [3] and 
[4]), and the close agreement between the Qjct — ^radio 
relations is puzzling, given such a large predicted offset. 

However, the radio luminosity is affected by a num- 
ber of other source parameters, including the density 
of the medium into wh ich the radio source expands 
(jBarthel fc ArnaudI fl99l . For a given radio luminos- 
ity, Qjet oc (Paper 2; Willott et al. 1999), where p 
is a characteristic ambient gas density. FRI radio sources 
typically inhabit more dense environments than FRII ra- 
dio sources (e. g. iZirbellllQQTt iMiller et al.lll999[ ). so the 
environment dependence will reduce, to some extent, the 
predicted offset in radiative efficiency described above. 

We test this hypothesis by comparing the position of 
Cygnus A in the Qjet— -^151 plane (Figure[3]) to that of our 
sample of FRII radio galaxies. Cygnus A is known to lie 
in a high density environment, more similar to the cluster 
environments of the FRI radio galaxies than the group 
or field environments of typical FRII radio galaxies, and 
therefore its luminosity is significantly "environmentally 
boosted" relat ive to similar sources loca ted in less dense 
environments (jBarthel fc ArnaudI 1X9960 . In Figure H it 
can be seen that Cygnus A is indeed an outlier compared 
to our sample of FRII radio galaxies, shifted to higher ra- 
dio luminosity. This offset is due to the "environmental 
boosting" of Cygnus A resulting from its dense environ- 
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ment. iBarthel fc ArnaudI ()1996D estimate that the himi- 
nosity of Cygnus A would be reduced by up to a factor of 
30 if it were located in a field environment. This would 
shift Cygnus A to align with the other FRII radio galax- 
ies in the Qj^t — ^radio plane. However, Cygnus A is 
only a factor a few below the FRI Qjot — ii-adio relation, 
and so the environmental dependence cannot be the only 
compensating factor. We note that the effect of "environ- 
mental boosting" in dense environments may be counter- 
acted to some degree by a possible positive correlation 
between the fraction of energy i n non-radiating parti- 
cles and the environment density (jHardcastle fc CrostonI 
[20T0t ICroston et al.l[20Tll) . 

Strong shocks driven by powerful FRII radio galaxies 
may also decrease their radiative efficiency relative to 
FRI radio galaxies. Unlike FRI radio galaxies, power- 
ful FRII radio galaxies drive strong shocks that sweep 
up and heat the am bient gas ahead of the cocoon (e. g. 
ICroston et aDl20Tll) . The energy associated with these 
shocks can be a signific ant, indeed dominan t, factor in 
t he FRII energy budg et (IWorrah et al.l[20T2h . 

IWillott et al] ()1999i ) argued that the Qjot — i^radio re- 
lation should be independent of source size. In Paper 2 
we show that this is in fact not the case, due to signifi- 
cant energy losses suffered by radiating electrons though 
inverse Compton scattering. These losses result in an ef- 
fective steeping of the spectral index a in Equation 1131 
and therefore a decrease in radiative efficiency. For the 
oldest sources the jet power can easily be underestimated 
by as much as a factor of three. 

Finally, we note that the X-ray cavit y jet power 
estimates used by iCavagnolo et al.l ()201C1[) . and hence 
the normalisation in Equation 1121 may be underesti- 
mated. Shocks, which are currently ignored in the 
X-ray cavity jet power calculations, may be energeti- 
cally important in FRI radio galaxies for a substantial 
part of their evolution, at least in some sources (e.g. 
McNamara ct al. 2005; Fab ian et al.ll2006l : iForman et al.l 
2007tiWise et al.ii2007c iBirzan et al .1120081 ). Furthermore, 
the buoyancy timescale used to estimate the source age 
could be an overestimate of the true source age, re- 
sulting in systemat ically underestimated jet power mea- 
surements (Wise ct "alll2007l: iMcNamara fc Nulsenll2007l: 
IBirzan et al.u200a) . 

A somewhat surprising corollary of our results, and 
the above discussion, is that an extrapolation of the 
Qjet — ^radio relation for low power FRI radio galax- 
ies provides a reasonable approximation for high power 
sources, despite their vastly different lobe energy bud- 
gets. 

5.3. The slope of the Qjct ~ iradio relation 

The exponent in Equations [TT] and [12] is < 0.7. In 
contrast, the predicted exponent in the Qjct — -^radio re- 
lation (Equation [T5|) . for a typical spectral index a w 0.8, 
is approximately 0.8. 

However, in a flux limited sample, selection bias will 
cause a systematic increase in radiative efficiency (a de- 
crease in the Willott / factor) with increasing radio lu- 



minosity: the low luminosity end will be populated by 
sources with both low kinetic power and low radiative 
efficiency (high /), while the high luminosity end will be 
populated by sources with both high kinetic power and 
high radiative efficiency (low /). For this reason alone, 
it is clear that a single value for the / factor cannot be 
applied to the entire source population, and the selec- 
tion bias will flatten (reduce) the observed slope of the 
Qjct — iradio relation, potentially accounting for the dif- 
ference between predicted and observed slopes. 

Relatedly, we note that given the common assumption 
of / = 20 based on the X-ray cavity jet power measure- 
ments for FRI radio galaxies, the Willott et al. relation 
(Equation [T]) predicts that the most luminous FRII radio 
sources (L151 ^ 3 x 10^* W Hz~^ sr~^) have jet power in 
the order of lO''^ erg s~^. This is equivalent to the Ed- 
dington luminosity of a 10^° solar mass black hole, and an 
order of magnitude greater than estimates of jet power 
in sa mples of radio galaxies (e.g. iRawlings &: SaundersI 
llQQlt) . It follows that / < 20 for the highest luminosity 
bins. 

6. CONCLUSIONS 

We have presented a new method to measure the ki- 
netic power in AGN jets based on the observed size and 
luminosity of the jet terminal hotspots. With this new 
method we were able to confront, from a new perspective, 
an emerging scenario in which the fraction of energy in 
non-radiating particles differs between the two morpho- 
logical clas ses of radio galaxy (F RI /FRII). IBirzan et al.l 
(2008) and ICroston etall (|2008D estimate that the ratio 
of energy in non-radiating particles to the energy in ra- 
diating particles, k, could be as high as 4000 in some 
sources, with a median value o f app roximately 200. In 
contrast, ICroston et al.l (|2004 120051 ) demonstrate that 
fc < 1 in FRII radio galaxy lobes. Such a large difference 
in the lobe energy budgets suggests a large difference be- 
tween the radiative efficiency of the two morphological 
classes, of more than an order of magnitude. To test this 
hypothesis, we estimated the jet kinetic power using the 
new method based on observed hotspot parameters, for a 
carefully selected sample of FRII radio galaxies. We com- 
pared the resulting Qjct — ^radi n relation to that deter- 
mined for FRI radio galaxies bv ICavagnolo et all ()2010f ) 
based on jet power measurements determined using the 
X-ray cavities method. 

We find approximate agreement between the FRI and 
FRII radio galaxy Qjet — iisi relations, which is os- 
tensibly difficult to reconcile with the differing lobe en- 
ergy budgets. However, a combination of environmen- 
tal factors, spectral ageing and strong shocks driven by 
powerful FRII radio galaxies reduces the radiative effi- 
ciency of these objects relative to FRIs, and conspires to 
move them onto the FRI Qjct — ^151 relation. An un- 
expected outcome of our work is that an extrapolation 
of the Qjct — iradio relation determined for low power 
FRI radio galaxies provides a reasonable approximation 
for high power sources, despite their apparently different 
lobe energy budgets. 
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APPENDIX 

THE VALIDITY OF THE ASSUMPTION Fe,jet ~ C x Fm.jet 

In this section we consider the requirements for the jet Lorentz factor in order to apply the assumption -pEjot 
c X Fivijct- Consider the ratio 



cx Fm 



jet 



E,jet 




(Al) 



Let us assume a tangled magnetic field (isotropic distribution in solid angle) in the jet. Then {B'j_ — Bj^) ~ {1/3)B^ 
and (-6^) = (2/3)i?^. In that case the above expression can be re-written as follows. 

= P + (A3) 

In general, provided pc^ < ee± : we have ^ ~ 2 — 4, regardless of the relative energy densities of leptons, protons and 

magnetic field. However, if pc^ » ee±, we have f >> 1, in which case ~ P- 

Below, we plot the ratio c x i^Mjet/^Ejet as a function of bulk Lorentz factor F using Equation IA21 for various 

combinations of the parameters , , . 

iMullin fc Hardcastld ()2009[ ) found that the radio emitting plasma in kpc-scale radio jets has a characteristic bulk 
Lorentz factor in the range F = 1.18 — 1.49. However, inverse Compton modelling of kpc-scale quasar X -ray jets 
suggests bulk Lorentz factors in the order of F 10 (see e. g. [Kataoka & Stawarz 2005). Mullin & Hardcast lc (200(1) 
argue that this discrepancy between jet speed estimates indicates the need for velocity structure across the jet (i.e. 
spine-sheath type of model). In that case, the average Lorentz factor of jet material will be F > 1.2 — 1.5, and 
potentially much greater. Figure [6] illustrates that, at worst, the value of c x Fm will underestimate the true jet power 
by up to a few tens of percent. However, the most likely scenario is that c x Fm will be within a few percent of the 
jet power. Note that the estimate will be an underestimate, so that this method represents a lower limit to the true 
jet power. 



1.00 




Jet Lorentz Factor r 



Fig. 6. — Plot of the ratio of speed of light times the momentum flux in the jet to the energy flux as a function of jet Lorentz factor for 
a range of jet plasma conditions. 
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